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Abstract. This course is an introduction to the physics of braneworlds. We concentrate on braneworlds with only one extra- 
dimension and discuss their gravity. We derive the gravitational equations on the brane from the bulk Einstein equation and 
explore some limits in which they reduce to 4-dimensional Einstein gravity. We indicate how cosmological perturbations from 
braneworlds are probably very different from usual cosmological perturbations and give some examples of the preliminary 
results in this active field of research. 

For completeness, we also present an introduction to 4-dimensional cosmological perturbation theory and, especially its 
application to the anisotropics of the cosmic microwave background. 



1. INTRODUCTION 

During recent years, cosmology has become one of the most successful fields in physics. The precise measurements 
of the anisotropics in the cosmic microwave background have confirmed a simple 'concordance model' : The Universe 
is spatially flat. Its energy density is dominated by vacuum energy (or a cosmological constant) which contributes 
about 70% to the expansion of the Universe, Q.\ ~ 0.7. The next important contribution is cold dark matter (CDM) 
with £lcDM — 0.3. The expansion velocity of the Universe is given by the Hubble constant Hq = 100/ikm/sMpc, with 
h ~ 0.7. Baryons only contribute a small portion of ~ 0.02. Massive neutrinos contribute similarly or less. 

The structures in the Universe (galaxies, clusters, voids and filaments) have formed out of small initial fluctuations 
which have been generated during inflation and have an almost scale-invariant spectrum, n = 1 ± 0. 1 . There is probably 
also a small amount of tensor fluctuations (gravity waves) generated during inflation, which however has not yet been 
detected. 

All the above numbers are accurate to a few percent and will be measured even more precisely with ongoing and 
planned experiments. This situation is unprecedented in cosmology. About twenty years ago, these numbers where 
known at best within a factor of two or even only by their order of magnitude. The concordance model is in agreement 
with most cosmological data, most notably the CMB anisotropy measurements, supernova type la distances (see 
contribution by Varun Shani), statistical analysis of the galaxy distribution, constraints from cosmic nucleosynthesis, 
cluster abundance and evolution etc. 

However, on a theoretical level our understanding has remained poor. We have no satisfactory answers to the 
questions: 

• What is dark matter ? 

• What is dark energy? 

• What is the 'inflaton'? Or what is, more precisely, the physics of inflation? 

• How can we resolve the Big Bang and other singularities of classical general relativity? 

There is justified hope that the last question could be resolved within a theory of quantum gravity, which is anyway 
needed if we want to put all fundamental interactions on a common footing. At present, the most successful attempt 
towards a theory of quantum gravity is string theory. This theory is based on the assumption that the 'fundamental 
objects' are not particles but one dimensional strings. Particles then manifest as excitations, proper modes, of strings. 
It would lead us much too far to give an introduction to string theory at this point. The interested student will have to 
study the two volumes of Polchinski [1]. 

In this course we next give an introduction to braneworlds, or, more generally, to physical effects of extra dimen- 
sions. In Section 3, we derive the gravitational equations for braneworlds with one co-dimension from Einstein's 
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equations in the bulk. We then discuss in detail the Randall-Sundram II model, its background and its perturbations. 
In Section 5, we give an introduction to 4-dimensional cosmological perturbation theory and, especially to the CMB 
anisotropy spectrum. Only after this we are ready for braneworld cosmology in Section 6. We write down the most 
general brane cosmology in an empty bulk and we investigate some of its modifications w.r.t. 4-dimensional cosmol- 
ogy. In particular, we discuss the modification of the slow roll parameters during braneworld inflation. We also present 
one example of the modifications in the evolution of cosmological perturbations which are relevant in braneworlds. 
We end with some conclusions. 

Notation: We use capital Latin indices A,B, - ■ to denote bulk coordinates, lower case Greek indices IJ.,V,- - for 
coordinates on a four dimensional brane, and lower case Latin indices, /, /, • • • for spatial 3-dimensional quantities. 
We sometimes also use bold symbols to denote 3d spatial vectors. We use the metric signature (—,+,••• , +). The 
4-dimensional Minkowski metric is denoted by (J7^v)- 

Throughout we set c = 7i = ^Boitzmann = 1 so that time and length scales are measured in inverse energies (usually 
GeV's), and mass and temperature correspond to an energy. The four dimensional Newton constant is then given by 
G4 = 0.67 X lO^^GeV 2. Useful relations in this set of units are 1 = 0.2GeV fm and 1 eV = 1.16 x lO'^K. Here fm 
= femtometer = 10^ ^^m. 



2. BASICS OF BRANEWORLDS 
2.1. What are Braneworlds? 

The interest of string theory Ues in the fact that it may provide a unified description of gauge interactions and 
gravity. Its weak point is, that it is extremely hard to make predictions from string theory which are testable at energies 
available in experiments. The reason for that is that string theory probably fully manifests itself only at very high 
energies of the order the Planck scale. The observed 4-dimensional Planck scale is given by Newton's constant, G4. In 
our units with h = c=l the Planck scale is E4 = M4 = 1 / v^47rG4 ~ 3 x 10'**GeV. This energy scale cannot be achieved 
by far at terrestrial accelerators (the LHC presently under construction at CERN will achieve about 7000GeV). 

Nevertheless, string theory makes some relatively firm predictions which might lead to observational consequences 
at low energy. First of all, it predicts that spacetime is ten-dimensional with one time and nine spatial dimensions. 
Since the observed world has only four dimensions, one usually assumes that the other six are compact and very small, 
so that they cannot be resolved by any physical experiment available to us so far. 

Furthermore, string theory predicts the existence of so called p-branes, p+ 1 -dimensional sub-manifolds of the 
ten dimensional spacetime on which open strings end. Gauge fields and gauge fermions which correspond to string 
end points can only move along these p-branes, while gravitons which are represented by closed strings (loops) can 
propagate in the full spacetime, the 'bulk' . 

This basic fact of string theory has led to the idea of braneworlds: it may be that our 3 -|- 1 -dimensional spacetime 
is such a 3-brane. If this is so, only gravity can probe the bulk and the additional dimensions can be much larger than 
the smallest length scale which we have probed so far, which is of the order of (200GeV)"' ~ 10" '^m. Actually, 
Newton's law has been tested only down to scales of about 0.1mm [2]. Hence, in the braneworld picture where only 
gravity can probe the extra-dimensions, these can be as large as 0.1mm= 10~^m. In the next subsection we show how 
this fact can be employed to address the hierarchy problem. 



2.2. Lowering the fundamental Planck scale 

The fact that the 4-dimensional Planck scale, M4 ^ lO'^GeV is so much larger than the fundamental scale in 
elementary particle physics, the electroweak scale Eew ~ lO^GeV is called the hierarchy problem. Apart from it 
seeming unnatural to have two so widely separated scales to describe fundamental physics, a more serious problem 
is the fact that as soon as we have a unified quantum theory which describes also gravity, the scale M4 will enter in 
quantum corrections of all electroweak scale quantities which are not especially protected e.g. by symmetries and it 
will therefore completely spoil the so successful low energy standard model. 

Here we show that within the braneworld picture, it is possible that the 4-dimensional Planck scale is not funda- 
mental but only an effective scale which can become much larger than the fundamental Planck scale Mp if the extra- 
dimensions are much large than Mp^. Our argument goes back to Arkani-Hamed , Dimopoulos and DvaU (1998) [3]. 
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Let Mp be the fundamental Planck scale and L the size of n extra dimensions. In addition there are 3 large spatial 
dimensions (and time). For simplicity we assume the n extra-dimensions to be rolled up as a cylinder, (5')®" iS^M?. 
The gravitational constants G(„_|_4) and G4 are defined by the force laws of gravity. Due to the Gauss constraint these 
must have the forms 

and F4 = G4 — 7j — . 



f-n+l 



A" 



7\ 



\J_ 



FIGURE 1. A 1 + 1 dimensional cylinder around a 1-dimensional grid of mass points, corresponding to the de-compactiflcation 
of a mass point on a circle. 

On small scales, r <C i, an observer on the brane sees n + 4 dimensional gravity, while on large scales, r^L, the 
cylinder can no longer be resolved and simple 4 dimensional gravity is observed. In order to relate the constants G4 
and G(„_|_4), we de-compactify the compact dimensions leading to an n-dimensional lattice of masses which looks from 
far like a hypersurface of mass density m/L" (see Fig. 1). Around the mass distribution, we now form a cylinder C of 
dimension n + 2, length d and radius r (see Fig. 1). To satisfy Gauss' law, we require 

Fj^da = Sf2+n)G(4+n) X (mass in C). 



'(2+«)"-^(4+«) 

Here, F± is the component of the force normal to the surface and Sj is the volume of a j dimensional sphere, 
Sj = 2n^j+^y^/r{[j + l]/2), where F denotes the Gamma-function, F([; + l]/2) = {[j - l]/2)! for integer values 
of U+l)/2. 

The first integral is 4nr^d"F{r) while the mass inside the cylinder is md"/L'^. With the 4-dimensional force law this 
implies 



G4 = 



'^(2+«) G(4+„) 



(1) 



4k L« 

In order to relate the gravitational constant to the Planck mass, we express Newtonian gravity in terms of an action 
principle. For a static weak gravitational field ^ and a mass density p in 4 -|- n dimensions, we can obtain the Poisson 
equation by varying the action 



M 



i2+n) 



(4+«), 



+ ... 



Integrating out <j), we obtain again the Newtonian force law and the relation 

r 1 

U2+n) _ '^(4+») 

^{2+n) 

Integrating the Lagrangian over the extra dimensions relates the 4- and 4 + n-dimensional Planck masses by 



(2) 



Together with Eq. (2) this reproduces the result (1). For a sufficiently large length scale L, M4 ~ 3 x lO^^GeV can 
therefore be much larger than the fundamental higher dimensional Planck scale Mp. 

Experimental 'micro-gravity' bounds [2] require L < 0.1mm. For n = 2 and L ~ 1mm the fundamental Planck scale 
can be of the order of the electroweak scale, Mp ^ (1 — 10)TeV . This Planck scale seems to be in agreement with 
most other bounds (cooling of supernovae, evolution of the Universe, etc) but leads to the very interesting prospective 
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that effects from string theory might be observable at the Large Hadron Collider (LHC) presently under construction 
atCERN [41. 

Using 'large' extra-dimensions, L ^ \ the fundamental Planck scale can therefore be of the same order as 
the electroweak scale. On the other hand, there is no explanation for the length scale L ~ 0.1mm~ 10~^eV. So the 
hierarchy problem has not actually been solved, but it has been moved from an energy hierarchy to an unexplained 
length scale. The hope is, that there would exist solutions of string theory which lead to such a scale dynamically. 



2.3. New Physics from higher dimensions 

Kaluza- Klein modes. We denote the four brane coordinates by and the additional n bulk coordinates by y". For 
simplicity, we consider a bulk where the extra dimensions are rolled up in a cyhnder of circumference L. In the general 
situation where the compact extra dimensions form a non-flat Calabi-Yau manifold, ^, the exponentials below have 
to be replaced by the corresponding eigen-functions of the Laplacian on ^. The case of a warped geometry, where the 
extra dimensions may even be non-compact, will be discussed separately in Section 4. 

Be now a massless scalar field in the 'bulk', (j>{x,y) with <!){■■■ ,y''+L---) =<j){--- , j", • • • ). We can expand the y 
dependence of in Fourier series 

(l>ix,y) ^Y.^jix)exp{i27tj-y/L) 
j 

Since satisfies the massless wave equation, V^^^^ = 0, in the bulk, for a 4d observer which cannot resolve the scale 
L, the modes j will become massive fields, 

(V^ + m])^j = with m] = . 
These modes of the gravitational potential give raise to exponential corrections to Newton's law. 

If the extra dimensions are large, the first few masses can be very low, but if the graviton is the only bulk field, its 
weak coupling to other bulk modes leaves the theory nevertheless viable [4]. 

As we shall see in Section 4, the mass spectrum can even be continuous, < m < °o, if there are non-compact extra 
dimensions. 



Higher dimensional spin modes. The spin states of massless particles in d space time dimensions are characterized 
by an irreducible representation of SO{d — 2). For d = A, massless particles always carry a 1-dimensional represen- 
tation of 50(2) Taking into account also parity, this leads to the two helicity modes of all massless particles in 
4 dimensions, independent of their spin. This situation changes drastically if we allow for extra dimensions. Let us 
consider = 5: a massless particle of spin smd = 5 spacetime dimensions, carries the representation of 50(3) 
and thus has 2* -h 1 spin states, like a massive particle of spin * in 4 dimensions. 

Let us, for example, consider the graviton. In 4 -|- 1 dimensions it has the 5 heUcity states of the tensor representation 
of SO(3). Projected onto a 3 + 1 brane, two of them become the usual spin 2 graviton, two are a spin 1 particle, a gravi- 
vector and one has spin 0, the gravi-scalar. 

The gravi-vector couples to the ;U4 components of the energy momentum tensor. Interpreting these as the electro- 
magnetic current and the gravi-vector as the electromagnetic potential A^, the five-dimensional Einstein equations 
lead to Maxwell's equations for and J^. This is the so called Kaluza-Klein miracle, which is also true if any, 
non-Abelian gauge group replaces the one-dimensional torus which plays here the role of the electrodynamic gauge 
group f / ( 1 ) . This finding of Kaluza and Klein [6] has evoked an interest in extra-dimensions in the 20ties, long before 
string theory. 



' The full little group for massless particles is actually IS0{2), the group of two dimensional Euclidean motions and rotations. But for finite 
dimensional representations of 750(2) the translations are acting trivially. The reason why not all representations of the universal covering group 
of SO(2), R, have to be considered, namely e'*'', £ R is rather subtle. On the classical level, where SO(2) is the relevant group, this problem 
disappears. A full discussion of the finite dimensional representations of the Poincar^ group for massless particles in d = 4 can be found e.g. in [5]. 
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The positive aspects of the vector sector are, however, over shaded by the problems coming from the gravi-scalar. 
It couples to the four-dimensional energy momentum tensor and modifies gravity. It leads to a scalar-tensor theory of 
gravity with several observable consequences. For example, one can calculate the modification in the slowing down 
of the binary pulsar [7] PSR1913H-16 due to the radiation of gravi-scalars. In Ref. [8] it is shown that in the simple 
case of a 5-dimensional cylindrical bulk, this leads to a modification of the quadrupole formula by about 20%, while 
observations agree with the quadrupole formula to better than 5%. If there are more than one extra-dimensions, there 
are several gravi-scalars and this problem is only enhanced. 

Clearly, a modification of higher dimensional gravity is necessary to address the problem. Usually, one gives a mass 
to the gravi-scalar. There are several ways to do this and the resulting four dimensional theory in general depends on 
this choice. One proposal is the Goldberger- Wise mechanism [9]. 

Another solution is offered by non-compact extra-dimensions. As we shall see, in curved spacetimes, extra- 
dimensions can even be infinite. Due to a so called 'warp factor' they become very small when seen from the brane. For 
infinite extra dimensions it can happen that the gravi-scalar represents a non-normalizable and therefore unphysical 
mode. This is precisely what happens in the RandaU-Sundrum model and we shall discuss it in this context in Section 4. 



3. GEOMETRY OF FIVE DIMENSIONAL BRANEWORLD GEOMETRY 

From now on we restrict ourselves to five dimensional braneworlds, i.e. braneworlds with only one extra-dimension. 
The idea here is still that spacetime has 10 (or for M-theory 11) dimensions, but 5 (or 6) of them are compactified 
to a static manifold of about Planck scale, while one of them is large. This picture is motivated mainly from 11- 
dimensional M-theory, e.g. the Horava-Witten model [10], but we shall not try to implement a realization of this model 
here. Nevertheless, it is important to note, that one co-dimension differs significantly from more than one. The main 
point is that the 3-brane spUts space into two parts, the 'left' and the 'right' hand side of the brane. We shall see, that in 
this case it is possible to determine the gravitational equations on the brane by simply postulating Einstein's equations 
in the bulk. This is no longer possible in the case of two or more extra-dimensions. 

In this section we derive and discuss these brane gravity equations. In the next section we shall apply them to the 
Randall-Sundrum model, which we consider as beeing so far the most promising braneworld model. 

To determine the gravitational equations on the brane, we start from the basic hypothesis that string theory predicts 
Einstein gravity in the bulk, 

Gab = K5TAB ■ (3) 

We want to discuss in detail the situation of a 3- brane in a 5-dimensional bulk. We denote the bulk coordinates 
by (x^) = {x^,y), where (x'') are coordinates along the brane and yis a transverse coordinate. We denote the brane 
position by y = yh, in general yt, depends on the point on the brane, yt, = ybix^)- A more general embedding of the 
brane will be discussed below. Very often we consider an energy momentum tensor of the form 

TabHx'') = {x\y)) = ^§AB + 5^5^T^y{{x^))5{y - y,) . (4) 

Here A5 is a bulk cosmological constant, T^y is the energy momentum tensor on the brane and K5 = 6n^G5 is the 
five-dimensional gravitational coupling constant. This is the most general ansatz if we do not allow for any matter 
fields in the bulk. 



3.1. The second fundamental form 

As above, gAB is the bulk metric. We denote the projection operator onto the brane by q^. The induced metric on 
the brane, also called the first fundamental form, is then given by 

SM^) =4{x^)q^y{x^)gAB{x^,yb{x^))- (5) 

We denote the co variant derivative in the bulk by Va and covariant derivative on the brane (with respect to the induced 
metric) by V^u. We also introduce the brane normal n, a vector field defined on the brane which is normal to all vectors 
parallel to the brane. Clearly, for an arbitrary vector field X = X^d^ along the brane, V^X ^ ^V^iX. The difference of 
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these two covariant derivatives is given by the extrinsic curvature also called the second fundamental form which we 
now introduce. 

Be X = X^d^ and Y = ¥^^3^ two vector fields on the brane. Their covariant derivative on the brane is given by 



WYX = {Y^'^^,X' + ^lpX^'YP)^y, 
while the covariant derivative in the bulk is 

"VyX = {Y^'^^X' +TlpXnP)dy + Tlpd4 
Therefore, there exists a bi-linear form K^y on the brane such that 

'"VYX = VYX + K{X,Y)n . 



(6) 



Since ''VyX - ''VxY = [Y,X] is tangent to the brane, K{XJ) - K(Y,X) = 0, hence K is symmetric. K is called the 
'second fundamental form' or 'extrinsic curvature' of the brane. Its sign is not uniquely defined in the hterature; we 
shall use Eq. (6) as its definition. 

Since n is the unit normal of the brane, g{n,^'VYX) = K{Y,X). But as n is normal to the brane vector field X we 
haveO = ''WY{g{n,X))=g{''VYn,X)+g{n,''WYX), so that 



K{Y,X) = -g{^VYn,X) 
In components, K{X,Y) = Kf^yX^Y^, this becomes 



g{''VYn,X)+g{''Vxn,Y) 



(7) 



Close to the brane we can choose coordinates such that 

gAsdx^dx^ = dsl = g^ydx^dx^ +dy^ . 

In these so called 'Gaussian normal coordinates', n = dy and ^Vj^rtv = r^v = ~(1 l'^)SnvA- The second fundamental 
form then becomes simply 

^iiv = —ijdyg^y . 

For general coordinates (z^ ) on the brane we have to define a brane parameterization x^ =X^{z^). The vector fields 
(e/i) = {d^X^{z)dA) then form a basis of tangent vectors on the brane. In terms of these one obtains by means of 
Eq. (7) 

Kny = -\ [gAB{epyrfi + 4d^rfi) + e^elrfgAB^ , (8) 



where a connma denotes an ordinary derivative, = 



3.2. The junction conditions 

Einstein's equations with a thin hyper-surface of matter become singular since there is a 5- function in the energy 
momentum tensor. Integrating them once across the brane leads to the so called junction or jump conditions (of Israel, 
Lancos, Darmois, Misner) [11, 12, 13, 14]. 

Before we come to the algebraically more complicated situation of general relativity, let us first recall the well 
known junction conditions of electrostatics: we consider a conducting boundary surface (e.g. capacitor plate) with 
surface charge density p and current density j along the surface. 

The homogeneous Maxwell equations require that the tangential part of the electric field and the normal compo- 
nent of the magnetic field, B±, are continuous across the boundary. This is usually derived by the following argument: 
denoting the two sides of the capacitor plate by the super-scripts -|- and — the homogeneous Maxwell equations imply 
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boundary 



FIGURE 2. A conducting boundary (bold line). We indicate the integration surface 5 normal to it and the 'pill box' type volume 
V encompassing both sides. Our equations become exact in the Umit where the height of the pill box and the width of the surfaces 
5 approach zero. 



for some surface S spanning from one side to the other of our boundary or for a volume V encompassing a Uttle of 
both sides of the boundary (see Fig. 2) 



= l^{VAE)da= JggE\\ds= L{E+-Ep hence 

[^ll] = E+-E^= 0, and 

= fy{VB)dv^ JgyBda= S{B+-Bl) hence 

[5J = B+-Bl= 0. 



(9) 



In the same manner, integrating the inhomogeneous Maxwell equations impUes for the normal component of E and 
the tangential component of B 

[Ej_]=El-E- = 4np, (10) 
[B||]=B^-B^ = AnjAn. (11) 

Similar junction conditions exist also for Einstein's equations. (Lanczos, 1922, Darmois 1927, Misner & Sharp 
1964, Israel 1966, see Refs [11, 12, 13, 14]), the so called junction conditions. To obtain them, we have to split the 
geometrical quantities into components parallel and transverse to a given hyper-surface. This is often also done in 4d 
gravity (3+1 formalism or ADM formalism), where one wants to study the time evolution of the metric on a 3d spatial 
hypersurface. Examples are numerical relativity, or canonical quantization of gravity where the canonical fields are 
the spatial metric components qij and their canonical momenta are given by the extrinsic curvature, Tiij = Kij. 

There one considers spacelike hypersurfaces, i.e. , hypersurfaces with timelike normal n, g{n,n) < 0. In braneworlds 
we have a timelike hypersurface with spacelike normals, g{n,n) > 0. 

Using a sUcing of spacetime into 4d hyper surfaces, one can express the 5d Riemarm curvature in terms of the 4d 
one and the extrinsic curvature. The equations are relatively simple if we write them in Gaussian coordinates (Gauss- 
Codazzi-Mainardi formulas, see e.g. [14]) 

Xa/3 = Kap+^yoi^p-^yP^« (Gauss formula) (12) 
^^\va = ^vK^a-'^aKnv (Codazzi formula) (13) 

^^m4v = dyKfiv+K^pK^ (Mainardi formula) (14) 

= 'R^'\y + 2'R'\^=R + 2dyK-K^-K„pK"P . (15) 
For the last equation we have used dyg^^ = 2K^^ so that 

g^'dyK^, = dyig^^'K^y) - {dygnK^., = dyK - IK^^K^y , K = xf, . 
From these we can determine the 5-dimensional Einstein tensor, 

'Gt = ^-[-R + K'-K^'^Ka^] (16) 
= Vf^K-VyK; (17) 



5g(J = G'i+2K^K^-KK^^ + ^yK^^-5l^^yK+^g^y{K^+K„pK"P). 



(18) 



The derivatives wrt y indicate that in order to determine the 5-dimensional Einstein tensor, it is not sufficient to 
know the second fundamental form on the brane itself, but we also have to know it on both sides of the brane. 
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We now derive equations on the brane from the bulk Einstein equation, Gab = ^sTab which we assume to be valid 
as a low energy consequence from string theory. To identify the energy momentum tensor on the brane which contains 
a delta-function in j-direction, we define 



Si = ]im 



/ Tidy 



^B 

'Jb 
so that 

lim / * ^(Ady = KsSi . 

The 4d metric is continuous and also K^v has no delta-function in y, but possibly a jump across the brane. This means 

c 4 c 4 c 

that G4 and have no delta-function. Only G^v may have one stemming from the term dyK/^y — gpLvdyK if K^y 
has a jump. Hence as consequence from the junction conditions, we obtain the following relations for Sab 

= 5^ (19) 
= 5^ and (20) 
KsSl = [Kl]-5;[K] or (21) 

[Kj,] = K5{Sl-\5;S) where S = S^^=Si. (22) 



3.3. Z2 symmetry 

In addition to Gaussian normal coordinates (which one can always choose, at least locally) we now assume Z2 
symmetry: the two sides of the brane are mirror images. For the metric components this implies 

giiv{x^ ,yb+y) = gnv{x^-yh-y) (23) 
gnA{x^,yb+y) = -gpLiix^ ,yb-y) (24) 
gAA{x^,yb+y) = gu{x^,yb-y) (25) 

The same symmetry is required for Tab and any other bulk tensor field. 

Under this condition we have ^+ = —K- so that the junction conditions (22) reduce to 

2< = K-5(5(f-i5v^5). (26) 

If the braneworld satisfies Z2 symmetry, the brane energy momentum tensor determines the second fundamental form. 

However, we now show that even with Z2 symmetry, knowing the brane energy momentum tensor is not sufficient to 
determine the brane Einstein tensor. To demonstrate this we first rewrite the 4d Einstein tensor, for a general coordinate 
system in terms of the 5d Riemann tensor and the extrinsic curvature: 

G^y = ^GABq^qv +^ RABn^n^g^v — K^Kav + KK^y 

+ lgnv{K^-Kal}K"l^)-Eny . (27) 

Here ^RABn^n^ corresponds to ^/?44 = ^R4^4 in Gaussian coordinates and E^y = ^RABCon^n^g^q^ corresponds to 
^/?4^4v Eq. (27) in Gaussian coordinates is a simple consequence of the expressions (12,13,14) for the Riemann 
tensor. In a general coordinate system it can be found in Ref. [15, 16] (careful, the sign for K^y is different there). In 
5 dimensions the Weyl tensor is given by 

2 1 

RaBCD = 2 \gA[cRD\B - gBlcRD]A) " g^A[C^Z)]B^ + CabCD ■ (28) 

Here [AB] indicates anti-symmetrization in the indices A and B, and Cabcd is the 5-dimensional Weyl tensor defined 
by Eq. (28). It is easy to verify that Cabcd is traceless and obeys the same symmetries as the Riemann tensor, Rabcd- 
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Inserting the 5-dimensional Einstein Eq. (3) for ^Gab and Eq. (28) in the expression containing the Ricci tensor 
Rab, as well as in E^y, we obtain the 4d brane gravity equation 



+ 



TAB^qy + g^iv (TABn^'n" - -T) 



where 



EyLv — CABCon'^n'^ ifi 



(29) 



(30) 



is the 'projection' of the Weyl tensor along the brane normal. The Codazzi equation (13) gives in addition 

KsTABn^ql = V^/: - VyKl . (31) 

Because of the last term in Eq. (29), it is not sufficient to know the bulk energy momentum tensor and initial 
conditions for n, gAB and Kab to solve the gravitational equations on the brane. In addition we need to know E^y, 
components of the bulk Weyl tensor on the brane. The Weyl tensor, which is the part of the curvature which can 
be non-vanishing even if the energy momentum tensor vanishes, contains information on bulk gravity waves. Bulk 
gravity waves can flow onto, respectively be emitted from the brane and thereby affect its evolution. This information 
is encoded only in the full bulk initial conditions. Therefore, to determine the evolution of the brane matter and 
geometry, in principle we have to solve the full bulk equations ! Only in situations with very special symmetries this 
can be avoided. However, as soon as we want to perturb such symmetric solutions we have to take into account all 
the bulk modes, and we do expect the solutions to differ significantly from the results of 4-dimensional perturbation 
theory. 



3.4. Brane gravity with an empty bulk 

We now exemphfy the effect of the bulk Weyl tensor in the case of an empty bulk. This will be the situation which 
we study for the rest of these lectures. We assume that the bulk is empty up to a simple cosmological constant A5. 



Tab = -^gAB + qA^qB^ S^vS{y - yb) ■ 

\5 



(32) 



where S/^y is the energy momentum tensor on the brane. It consists of a brane tension A and the matter energy 
momentum tensor x. 



^lv■ 

The junction conditions read 

[Suv] 

Z2 symmetry requires that 



Siiv — ^gixv + "^v • 



first junction condition, 



second junction condition. 



3^Mv5 ) . 



(33) 

(34) 
(35) 

(36) 



Inserting our ansatz for Tab in the brane gravity equation (29) and using the second junction condition to eliminate 
the second fundamental form, we obtain 



with 



A4 = 
fCt = 



KfA/6 = 2/M| 
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and, as before denotes the projected 5d Weyl tensor, evaluated on either side of the brane (but not exactly on the 
brane where it may be ill-defined). The quantities K4 and M4 denote the 4-dimensional gravitational coupling constant 
and Planck mass respectively and T = T^^ is the trace of the matter energy momentum tensor. The relation between 
the 4- and 5-dimensional Planck mass in the braneworld approach is now obtained as follows: using ~ M^^ and 
K4 ^ M^'^, Eq. (39) shows that M| ~ M|/A. Using that the 4-dimensional cosmological constant is small, A4 IA5I, 
we have X ~ V-SAs/fCs ~ y/^ET^M^, so that M| ~ M^^-As = M^L with L ~ ^/-A5. 

In the limit Xx^v ^ f^af" we recover the 4-dimensional Einstein equation if E^v is negligible. The existence 
of this hmit depends crucially on the existence of a 4-dimensional brane tension. In order for the 4d gravitational 
coupling constant to be positive, the brane tension must be positive, X > 0. At high energy densities (in the early 
universe) the quadratic term a^y can become dominant and modify the dynamics (the expansion law of the universe). 
In general, there is an additional part, E^iv , carrying information from the bulk geometry and evolution, which can 
affect the brane evolution in a crucial way. 



3.5. Energy momentum conservation 

The Codazzi equation (13) together with Z2 synometry impUes 

V^,K-VyK; = K5T^ = 0. (41) 
With the Gauss equation (12) and Z2 symmetry this ensures energy and momentum conservation on the brane, 

VvT;=0. (42) 
From the 4-dimensional contracted Bianchi identities we obtain in addition 

VvZ?; = K-|Vvo-; . (43) 

Hence, the longitudinal part of E^v is fully determined by the matter content of the brane, while the transverse traceless 
part is not specified: 

£;,v=4v^^+4v (44) 



where WyEj^ ^^^^ = and 



4v = \ (Vm^v + Vv0m) with V^0^ = . 
Inserting this in Eq. (43), we obtain 



(45) 



For given initial conditions, this equation has always a unique solution 8^ on the brane which determines £'ja''- 

(TT) 

However, the transverse part, E^y is not determined by the brane energy momentum tensor; it comes from bulk 

(TT) 

gravity waves. Only if is^y = does the brane energy momentum tensor determine the brane Einstein tensor. As we 
shall see, already in quite simple situations this is not the case. 



4. THE RANDALL SUNDRUM MODEL 

We now consider an Anti-de Sitter (AdS) bulk, A5 < and would like to obtain Minkowski space on the brane. Since 
AdS is conformally flat, Ej^y = 0. A Minkowski brane can be achieved by setting T^jy = 0. If in addition the brane 
tension is related to the 5-dimensional coupling constant and the cosmological constant by 

X^kI/6 = -A5 , (46) 

Eq. (37) imphes G^y = 0. Eq. (46) is it the Randall-Sundrum (RS) fine tuning condition [17, 18]. Small deviations 
from the RS condition lead to an exponentially expanding/contracting brane. The 4-dimensional gravitational constant 
becomes 

A5 , A5 /-A5 

K4 = — -^ = -— >0, or, equivalently A = , and K4 = K5'i/— — . (47) 

6 A fQ V 6 
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4.1. The metric 



The following coordinates for (a part of) Anti-de Sitter will be useful for us: 

ds^ — e^^^'^l^r\^vdx^dx^+dz?' Gaussian coordinates (48) 
ds^ — [-] (7]^ydx^dx^ +dy^) \y\> conformal coordinates. (49) 



Einstein's equations. Gab — —^sgAB, give — The RS fine tuning requires X = y / . 

In their first model [17] (RSI model) Randall and Sundrum propose two branes, the first positioned at z = is 
called the hidden brane, and the second positioned at z — k£ is called the visible brane and represents our Universe. 
The gravitational force on the visible brane is suppressed by the factor exp(— 2A:) w.r.t. the hidden brane, leading 
to an enhancement by a factor exp(fe) of the apparent Planck mass. However, also this two brane model contains a 
gravi-scalar (also called radion) which has to obtain a mass by some non-gravitational mechanism. 




FIGURE 3. The RSI model with two branes. The visible brane (our Universe) a.tz = k£ and the hidden brane at z = 0. 




FIGURE 4. The RS2 model with only one brane at z = 0. 

This problem is resolved in the second model [18] (RS2 model). There, our universe is located on the brane at z = 
corresponding to y — £ and no second brane is present. The apparent 4-dimensional Planck mass as measured on the 
brane on scales much larger than £ is then again given by 



K4_ 

As 



3M| 

£2 



:£-^K^ 



Furthermore, 



so that 



(50) 
(51) 
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Hence in the RS2 model, the AdS curvature scale I enters in the same way as L for cylindric Kaluza-Klein models. 
Since the scale I is limited by present day micro gravity experiments which have not detected any deviation from 
Newton's law [2], we have 

I < 0.1 mm implying A > (1 rey)"^, hence Mj = (m|/^)'/^ > 10^ JeV . 
4.2. Gravity waves in the RS2 model 

In order to see that the radion mode is absent in the non-compact RS2 model, we consider perturbations to the AdS 
metric. It is easy to show that one can always choose a gauge (local coordinate system) so that the perturbed metric is 
of the form 

ds^ = ^0 [-{\+m')dt^ -AZidtdx' + {{\-2^)5ij + 2Hij)dx'dx^ + AZidydx' (52) 

-A^dtdy + (1 + 2'£)dy^] (53) 

Here Htj and E,, S, are transverse (i.e. divergence free) and Hij is traceless. In other words, (5,E, = (9,S, = diHij = 
HI — 0. As we shall see there are 5 homogeneous modes in these 10 physical perturbation variables corresponding to 
the 5 gravity wave modes in 5 dimensions. We consider the homogeneous 5d wave equation in the bulk. 

Since the 3-brane is homogeneous, scalar-, vector- and tensor degrees of freedom decouple and we can consider 
them in turn (for more details see Section 5). 



4.2.1. Tensor perturbations 
We first consider only the tensor Hjj, so that the perturbed metric is given by 



ds'^=(^ {-dt^ + {5ij + 2Hij)dx^dx' +d/) . (54) 

We Fourier transform //,, in the 3-dimensional x coordinates and consider one mode with fixed wave vector k, so that 
Hij (f , y, x) — Hij (f , >') exp (;k • x) . Since spacetime is isotropic and homogeneous in x, different k-modes do not couple. 
The bulk Einstein equations, SGab = —^SgAs, for the Fourier mode k then give 

d^+l^- + Hij = . (55) 

The general solution to this equation is of the form 

Hij^h,neij with/i„ = e'""{myf[AJ2(my)+BY2{my)]. (56) 

Here is the polarization tensor, k'etj = e'^ = and (O^ = m^ +k^. The separation constant is arbitrary and 
can, in principle also be negative. J2 and Y2 are the Bessel functions of order 2. They are oscillating and decaying. 
Bessel functions represent "5-function normalizable" perturbations like harmonic waves in flat space, in the sense 
that [19, 20] 

h,nKn' — (wi ^ • (57) 

These are just the ordinary gravity modes of 4-dimensional mass m without a mass gap which are discussed in the 
original RS paper [18]. To find the correct weight l/>'^, we use that hm satisfies 



3 

y 

and thus h = h^/y^^^ satisfies the equation of motion of a scalar field in a flat 5-diniensional spacetime (with y- 
dependent mass term), 

15 

4y' 



a^ + d^--dy]h„ = 0, (58) 



°4 + ^'-£2 )'' = 0- (59) 
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This mode has to be normalizable w.r.t the Minkowski metric (no additional weight). 

As we have mentioned above, is arbitrary and can also be chosen negative. However, if < and therefore m 
is imaginary, it is more useful to decompose the two independent solutions in the form 

h„, = e'"'(|mb)2 [CK2{\m\y)+Dh{\m\y)] , (60) 

where K2 and I2 are the modified Bessel functions of order 2. Considering the behavior of the Bessel functions, one 
sees that I2 grows exponentially (see Fig. 5) and is clearly not normalizable {i.e. not square integrable with some 
weight which is a power law in y). Therefore, this mode is unphysical and we have to set D = 0. It is important to note 
that aP' =k'^ + can become negative in this case leading to (O = ±i\(o\. In other words, negative mass solutions 
become exponentially growing 'tachyonic' instabilities! It is still unclear whether these tachyonic modes are relevant 
for cosmological braneworlds. 




X 



FIGURE 5. The Bessel functions of order 2. 
In the limit m ^ only the l2-mode survives and we obtain /!(,„=o) = Cexp(a)f ) independent of the coordinate y. 
This zero-mode is normalizable with respect to the measure dy/y^. 
The general solution for a tensor perturbation is of the form 



h — ho 



hindm 



(61) 



At the brane position, y = yb = £, the perturbations must satisfy the junction condition (36). These represent boundary 
conditions for the perturbations Hjj in the bulk. On the right hand side of Eq. (36), we can in principle have an arbitrary 
perturbation of the matter energy momentum tensor. However, the only non-vanishing term of the tensor contribution 



to Tjjv is the traceless part of T/y, i.e. the anisotropic stress on the brane, 11,^ 

'2 



(T) 



5Ki, 



A short computation shows 
hence 



Vb 



2{dyH,j)\ 



where H*^^' are tensor-type anisotropic stresses on the brane. 

Let us first consider the homogeneous case II'^' = 0. For > 0, the solutions are of the form 

h = exp{±im){my)^ [AJ2imy) +BY2{my)] . 

The junction condition (62) then requires 



B = -A 



J\ {m£) 



'-{mlfA , 
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where the last expression is a good approximation for mi<^l. This is precisely the result of Randall and Sundrum [18]. 
It is not modified even if we allow for the negative mass modes, — > 0, because a physical solution has to be of the 
form 

h = Cexp{±tV-m^-k^){\m\yfK2{\m\y) , with dyh = \m\Cexp{±tV-m^-k^){\m\yfKi{\m\y) , (65) 

and since Ki has no zero, the junction condition (62) requires C = 0. 

But in a realistic brane universe, n^^) is not exactly zero. In cosmology, it is typically just a factor 10 smaller than 
other perturbations of the energy momentum tensor on the brane. We therefore cannot require C = 0. However, as 
long as n^^^ remains small, we do not expect the unstable modes to be present, so that C{k,m) = for < —m^. 
Within the framework of first order perturbation theory, the II^^' modes satisfy a Minkowski equation of motion and 
therefore they do not grow exponentially. Hence in this case, the A'-mode can only be excited for (O-^ = + np' > 0, 
i.e. > — OT^. However, it is not clear whether this remains true to second order, where the evolution of H feeds back 
in the equation of motion for H^^) . Actually, within a toy model, it has been shown that the full, non-linear evolution 
can be exponentially unstable even if the linear equations do not excite the unstable mode [21]. 



4.2.2. Vector perturbations 

We now consider vector perturbations only, so that, in generalized longitudinal gauge, the metric takes the form 



ds^ = {-dt^-4I.idtdx' + 5ijdx'dx^ +4Zidydx'+dy'^) 



(66) 



The bulk Einstein equations for a mode k of the vector perturbations Z and S are 

3 



d^--dy]Z={d^+k^)Z, (67) 



y 

^y-'-dy + ^)- = id? + '^)^, (68) 
dy-^^Z = -dtZ, (69) 

where Z and S are transverse vectors, (k • E) = (k • S) = 0. The constraint equation (69) fixes the relative amplitudes 
of E and S, showing that there is only one independent vector perturbation in the bulk (the "gravi-photon"). One can 
check that these equations are consistent, e.g. with the master function approach of Ref. [22]. 

As in the tensor case, the solutions are Bessel functions of order two (and one). Considering just one component 
Z = Z, one obtains the expected oscillatory modes for positive mass-squared, > 0, 



L = sxp{±i(ot){myf[AJ2{my)+BY2{my)] , (70) 
exp{±ia)t){myf[AJi{my)+BYi{my)] , (71) 



where fi) = \/rrP--\-k?-. These solutions have been found in Ref. [23]. For a negative mass-square, < 0, we obtain 
again tachyonic solutions. Like in the tensor case, the solution containing the modified Bessel function ly cannot be 
accepted as it is exponentially growing and thus represents a non-normalizable mode. However, the /Ty-solution is 
exponentially decaying and perfectly acceptable. For tachyonic vector perturbations with OJ^ = m'^ -|- A:^ < we have 

Z = C^xp(±\(ii\t){\m\yfK2{\m\y) , (72) 

S = ^Cexp{±\co\t)i\m\yfKi{\m\y) . (73) 

\m\ 

For large enough scales, —np- > k^, these solutions again grow exponentially. 
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The boundary conditions at the brane relate these perturbations to the brane energy momentum tensor. For the 
energy momentum tensor on the brane, the vector degrees of freedom are defined according to 



M=(l. ^mj-^M > (74) 

where V; and n -^^ are divergence- free vector fields and n -J^ = [i^^n -^^ + diYl^J'^] . The first junction condition simply 
requires that Z be continuous at the brane, which it is since the (modified) Bessel functions of even index are even 
functions. The second junction condition results in (for a detailed derivation, see [24]) 

dtE + dyL=K^V, (75) 
S = jqn^'') , (76) 
dtV = -A^n(^) . (77) 

The last equation follows from (75) and (76) and the bulk equations (67)-(69). It represents momentum conservation 
on the brane, which is guaranteed as long as we have vanishing energy flux off the brane and Z2-symmetry. 

Like for tensor perturbations, we consider homogeneous solutions, setting II^^) = V = 0. This requires S(|m|£) =0, 
hence 

B = -AprT^ form2>0, (78) 
Yi [ml) 

C=0 forw2<0. (79) 

Equation (75) is then identically satisfied. 

However, it seems more realistic to allow a small but non-vanishing anisotropic stress contribution n^^) and 
corresponding vorticity V. In this case, again, we can have solutions with C 7^ which can grow exponentially in 
time; hence small initial data can lead to an exponential instability like for tensor perturbations. 

Using the normalization condition (57) for the m — mode of the variable E °^ y (this is the one which enters as 
dynamical variable in the perturbed action, see [25]), one finds that / \Z\^/y^dy diverges logarithmically. Contrary to 
the tensor case, the vector zero-mode is not normalizable. Therefore, on the brane there is only the ordinary massless 
spin-2 graviton, but there are a continuous infinity of massive spin-2 and spin-1 particles (the modes discussed here, 
with m ^ 0). 

4.2.3. Scalar perturbations 

We now discuss the most cumbersome, the scalar sector. Scalar-type metric perturbations in the bulk are of the form 

£2 

ds^ = — \-{\+ 2V)dt^ - Amtdy 

y 

+{\-l(b)5ijdx'dx^ + {l+l'^)dy^] . (80) 
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The bulk Einstein perturbation equations for the mode k become, after some manipulations and introducing the 
combination T = ^ + »P (see Ref . [21]), 

^-^ = '^, (81) 

d^-hy^r={d^ + k^)T, (82) 

dy^+(dy^-]'^=-d,^ , (84) 



^(^dy-^^'^ = 3d^^ + k^{'^ + '^) , (85) 
3dt(dy^-^'rf^ =k^:^ , (86) 

dt {20-'^)= (^dy-^^^ . (87) 

Clearly these equations are not all independent, Eqs. (86)-(87) are identically satisfied if Eqs. (81)-(85) are. 
The so lutions a re obtained as for tensor and vector perturbations. For a positive mass-square, > 0, we find 
(© = \/m2+F) 

r = exp{±im){myf [A'j2{my) + B'Y2{my)] , (88) 
'^ = exp{±ia)t){myf[AJo{my)+BYo{my)] , (89) 

^ = ^exp{±i(Ot){myf [A'J2{my) +B'Y2{my) 

+AJo{my)+BYo{my)] , (90) 

»p = ^exp{±i(Ot){myf [A'j2{my) +B'Y2{my) 

-AJo{my)-BYo{my)] , (91) 

^ ±im^y^ . , . 

^ = — exp(±ia)f ) X 

2(0 

[{A' -3A)J,{my) + {B' -3B)Y,{my)\ , (92) 

2 2 

in tfi 
with A' = 3A „ ^ - , and S' = 3B , ^ , . (93) 

For a negative mass-square, < 0, we obtain (to = V — — k'^) 

Y = QX^{±m){\m\yfC'K2{\m\y) , (94) 

^ = tx^{±m)(\m\yfCKo{\m\y) , (95) 

1 , 
<I>= -exp(±(Bf)(|m|3') X 



[C'K2{\m\y)+CKQ{\m\y)] , (96) 

1 -J 
W = -exp(±cof)(|m|y)^x 

[C'K2{\m\y)-CKo{\m\y)] , (97) 

^= -^—^exp(±cot)[C' + 3C]Ki(\m\y) , (98) 
2a) 

I |2 

with C' = -3C , ,' ' , , (99) 
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where we have akeady used that the 7-mode is not normalizable and therefore cannot contribute. Like for vector and 
tensor perturbations, we find again tachyonic solutions with np' <0 which represent an exponential instability for 
sufficiently small wave numbers k (large scales). 

Determining the boundary conditions via the first and second junction conditions now requires a bit more care. Since 
we have already fully specified our coordinate system by the adopted choice of perturbation variables, we must allow 
for brane bending. We cannot fix the brane at yb = i, but we must allow foIy^^ =£ + S' and = —£— respectively. 
Fortunately, <^ is a scalar quantity and brane bending therefore does not affect vector and tensor perturbations. The 
anti-symmetry = —y^ is an expression of Z2-symmetry. The introduction of the new perturbation variable S{x^) 
describing brane bending enters the exressions for the first and second fundamental forms. From Eq. (5), we obtain 
q^y = g^v to first order, which implies that <i> and "P, hence have to be continuous. At the brane position, the 
perturbed components of the extrinsic curvature (8) are 



SKoo = \ 



S' 

0-3'P + 2- 



1 



5^0 = 



'P-3<I>-2- 



5ij + didj^. 



For the energy momentum tensor on the brane, we parameterize the four degrees of freedom according to 



(^Mv) 



(100) 

(101) 
(102) 

(103) 



where v,- = diV and n|J' = (didj — jASy) n^'^). With Eqs. (100)-(102), the second junction condition reads 



^ (2p + 3p) = <^-W+edt {dtS - 2^)+LdyW , 



3 

-An(^) 



= w-<^+edy<^ . 



(104) 
(105) 
(106) 
(107) 



Combining the time derivative of Eq. (105) with Eqs. (104), (84) and (107), we obtain momentum conservation on the 
brane, 

diV = ^An(^) + p . (108) 

Similar manipulations imply energy conservation on the brane, 

d,p = Av . (109) 

Like for tensor and vector perturbations, we look for solutions with vanishing brane matter. Setting 11^^^ = p = P = 
V = forbids brane bending, <^ = 0. Then Eq. (105) implies M{m£) = 0, thus 



B' -3B=-{A' -3A) 
C' + 3C = 



Ji{mt) 
Yi{mi) 
for OT^ < . 



for >0 



The other equations are all satisfied if we require separately 

B_B^_ J\{mt) 
A~A'~~Yi{mi) 
C = C' = for OT^ < . 



for m > , 



(110) 
(111) 

(112) 
(113) 
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Since B/A = B' /A', equations (1 10) are (1 12) are equivalent. 

As for vector perturbations, the m = scalar mode is not normalizable. Like for tensor and vector perturbations, we 
have found "scalar gravitons" which appear on the brane as massive particles. If the brane matter is unperturbed, only 
oscillating > solutions are possible. However, if we allow for non-vanishing matter perturbations on the brane, 
we can have C ^0 and the tachyonic modes < can appear exactly like in the tensor and vector sectors. 

It is not surprising that the same instabiUty appears in the scalar, vector and tensor sectors, because all modes 
describe the same bulk particle, the five-dimensional graviton. 



4.3. Green's function, correction to the Newtonian potential 



We want to determine the modification to Newton's law in the RS2 model. Since the extra dimension is not compact 
and there are massive (homogeneous) modes of all masses > 0, we expect a modification which is not exponentially 
suppressed. 

The 5d Green's function is defined by 

W^G{x,x) = 5^{x-x) , 

where is the 5d d'Alembertian in AdS spacetime and we have to glue together a Z2-symmetric solution on both 
sides which satisfies the homogeneous junction condition. We can obtain the retarded Green's function in the standard 
way from the homogeneous solutions of the equation (see e.g. [26]): 



+ 



dm- 



The first term comes from the m = solution and the functions Um are the properly normalized massive modes, 

fml Ji {mi)Y2 {my) - Yi {ml)J2 {my) 

The general retarded solution for a given energy momentum tensor on the brane is now of the form 

hfiv{x) = -2k5 j d'^x Gr{x,x')S^v{x') ■ 

For a stationary matter distribution it is simpler to use the Green's function of the spatial Laplacian which is related to 
Gr via integration over time 



G{x,y,x',y') 



dt' Gr{x,x') 



d^k 
(2^ 



,*-(x-x') 



k2 



■I 



dm 



um{y)um{y') 



OT2+k2 



4nr 2n Jo 



dmum{y)um{y') exp(-OTr) 



(114) 



(115) 



On the brane, y = y' = the first term gives the usual 1 /r behavior, r = |x — x' | . Expanding the second term to lowest 
order in i/r we obtain 

G{x,e,x',i) ~ ~^ 



4Kir 



1 



2r2 



This determines the Newtonian potential of a point mass on the brane with mass M, 



KiMGi 



-K4M 



1 



2r2 



(116) 



Since the extra dimension is non-compact, the correction is not exponentially suppressed but only as a power law. 
Away from the wall, the potential at large separation is given by 



2r2-h3/ 



87r(£-Fy)2 (r2+j2)3/2 



(117) 



The equipotential lines are shown in Fig 6. This formulas have been derived in Ref. [27] from which also Fig. 6 is 
drawn. 
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FIGURE 6. The equi-potential lines for the gravitational potential from a point mass on the brane in RS2. Figure from [27]. 

5. COSMOLOGICAL PERTURBATION THEORY IN 4 DIMENSIONS 

Before studying brane cosmology and 5d effects on cosmological perturbations, I present a brief introduction to 4d 
cosmological perturbation theory and some aspects of 4d cosmology. Much more details can be found e.g. in [28, 29]. 
Some knowledge of 4d cosmology is however assumed (Friedmann equations etc. as they can be found in the first 
chapter of standard textbooks on cosmology like Refs. [30] or [29]). Students who are familiar with this subject may 
skip this section. 



5.1. Perturbation variables 

The observed Universe is not perfectly homogeneous and isotropic. Matter is arranged in galaxies and clusters of 
galaxies and there are large voids in the distribution of galaxies. Let us assume, however, that these inhomogeneities 
lead only to small variations of the geometry which we shall treat in first order perturbation theory. For this we define 
the perturbed geometry by 

g^iv = gfiv + a^hiiv , gfivdx^dx^ ^ [—drf' + 'Yijdx'dx^) . (118) 

Here g^^ is the unperturbed Friedmann metric, a{r\) is the scale factor, Tj denotes conformal time and jij is the 
3d metric for a space of constant curvature K. The perturbations are assumed to be small, \h^y \ <^ 1. The energy 
momentum tensor is given by 

T>^ = T>i + e>;, Tl = --p, T)=P^ |0v^|/p«l. (119) 
The background energy density p and pressure p satisfy the Friedmann equations, 

JP^=(^\\k = ^fl2p + ifl2^4 (120) 

\a J 3 4 

p = -3(p + /5)jr , (121) 

where an over-dot denotes the derivative w.rt. conformal time rj. 

Without loss of generality we can choose the so-called longitudinal gauge so that perturbations of the metric are of 
the form 

^''^'^=\ B, (l-2<I>)^, + //,,j ■ ^122) 

Here B, is an divergence free vector and is a trace-free, divergence free tensor field. The scalar quantities ^ and <I> 
are called the Bardeen potentials. In the Newtonian approximation they are both equal and reduce to the Neewtonian 
potential. 
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We also decompose the perturbations into different 'Fourier modes', 
'P(j7,x)=Fk(xmJ?,k), <I.(j7,x)=Fk(x)0(j7,k), fi'(j7,x) = F^^)'(x)5(j7,k), Hij{r],x) =Y^'',^j{x)H{r],k) . 

In a Friedmann Universe with vanishing curvature, these are just ordinary Fourier modes, while in the general case, 
the functions are eigenfunctions of the spatial Laplacian with eigenvalue —k^. The functions Y^^' and F^^j 
correspondingly are vector- and tensor-type eigenfunctions of the spatial Laplacian with vanishing divergence. For 
later use we also define a scalar type vector and tensor as well as a vector-type tensor, 

^ -^-iViFk, y^^A=^-2(V,V;-l5i,A)7k, (123) 



€1 = ^(^^y^h^jy^)- (124) 



Av) 



Note that contrary to the vector-type vector field 7,- ', the vector field y,- ' is not divergence free. The same is true for 

the tensor fields Y^f and Ylf\ 

Let Ty =Ty+ 0y be the full energy momentum tensor. We define its energy density p and its energy flux 4- vector 
u as the time- like eigenvalue and eigenvector of Ty : 

T^u' = -pu^, u^ = -l. (125) 

We then parameterize their perturbations by 

p=p(l + 5), u = u^dt + u'di. (126) 

m" is fixed by the normahzation condition, 

MO = i(l-ip). (127) 

We further set 

u' = -v' = - ( VY^'^' + y(^)7(^)'") . (128) 

Here 5 is called the density contrast and (v') is the pecuhar velocity. 

We define Py = u^Uy -|- Sy , the projection tensor onto the part of tangent space normal to u and the stress tensor 

x^'" =P^P^T"P. (129) 

In the unperturbed case we have Tq = 0, Tj = p5j. Including perturbations, to first order we stiU obtain 

tO = tP = t^ = 0. (130) 
But Tj contains in general perturbations. We set 

T) = p[{l + nL)5} + n)], with n| = 0. (131) 
We decompose 11^ into scalar- vector- and tensor-type contributions. 

Another important variable is 

r = 7tL-—5 (133) 
w 

where = p/p is the adiabatic sound speed and w = p/p is the enthalpy. One can show that F is proportional to the 
divergence of the entropy flux of the perturbations. Adiabatic perturbations are characterized by F = 0. 
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We shall use also other perturbation variables describing the density contrast and peculiar velocity, which actually 
correspond to these perturbations in different coordinate systems (gauges). One can show that on sub-horizon scales, 
k ^ J^, on which perturbations are actually measurable, they all coincide. 

D ^ 5 + 3(l+w)(^)^, (134) 
Dg = 5-3(1 +w)4>, (135) 

a = y(^)-B(^), (136) 

n-yW = -z?W = a(^). (137) 

Here we use the customary name, a^^^ = —B^^\ for the vector-type metric perturbation. These variables can be 
interpreted nicely in terms of gradients of the energy density and the shear and vorticity of the velocity field [31]. 



5.2. Einstein's equations 

We do not derive the first order perturbations of Einstein's equations. This can be done by different methods, for 
example with Mathematica. We just write down the results. 



5.2. 7. Constraint equations 



AnGa^pD = -{k^-3K)^ (00) 
4nGa\p+p)V = k{jf^ + ^) (Oi) 

SKGa^{p+p)Q. = ^{2K-k^)G'^^'> (Oi) 



(scalar) 
(vector) 



(138) 
(139) 



5.2.2. Dynamical equations 



r(a)-'P) 

+ 



= SnGa^pU^^^ {i j) (scalar) 



k(a^^)+2(^]a^^^ 



(r) 



4KGa^ \^c^Dg + wr- -wn 
SnGa^pU^^'^ (ij) 
SnGa^pYl^'^^ (ij) 



(140) 
(141) 



{i i) (scalar) (142) 
(vector) (143) 
(tensor) (144) 



For perfect fluids, where IV- = 0, we have <t> = a^^'' <^ and H^'^^ obeys a damped wave equation. The damping 

term can be neglected on small scales (over short time periods) when 77^^ ^ 2K + k^, so that Z/*^^-' represents a 
propagating gravitational wave. For vanishing curvature, the scales kr\ \ are simply the sub-horizon scales. For 
A" < 0, waves oscillate with a somewhat smaller frequency, (0 = y/2K + k'^, while for A' > the frequency is somewhat 
higher than k. 
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5.2.3. Energy momentum conservation 



The conservation equations, VyT''^ = T-^ =0 lead to flie following perturbation equations. 

b^ + 3{c]-w){i)D, + {\+w)kV + 3w{i^T=Q 

^+(f)(l-3c?)^ = M^ + 3c?*) + ^l^. \ (scalar), (145) 



wk 



r-l i-f in 



a+(l-34(?)a = ^(.-f)nC) (.ec,„r)^ 



(146) 



These can of course also be obtained from the Einstein equations since they are equivalent to the contracted Bianchi 
identities. 

For scalar perturbations we have 4 independent equations and 6 variables. For vector perturbations we have 2 
equations and 3 variables, while for tensor perturbations we have 1 equation and 2 variables. To close the system we 
must add matter equations. The simplest prescription is to set Y = n,y = 0. These matter equations, which describe 
adiabatic perturbations of a perfect fluid give us exactly two additional equations for scalar perturbations and one each 
for vector and tensor perturbations. 

Another example is a universe with matter content given by a scalar field. We shall discuss this case in the next 
section. More complicated examples are those of several interacting particle species of which some have to be 
described by a Boltzmann equation. This is the actual universe at late times, say z ~ 10^". 



5.2.4. A special case 

Here we rewrite the scalar perturbation equations for a simple but important special case. We consider adiabatic 
perturbations of a perfect fluid. In this case n = and F = 0. Eq. (140) implies <t> = ^. Using the first equation of 
(138) and Eqs. (135,134) to replace Dg in the second of Eqs. (145) by and V, finaUy replacing V by (138) one can 
derive a second order equation for which is, the only dynamical degree of freedom 

* + 3 + c^)* + [(1 + 3c2) -K)-{l+3w) + + c^fe^Jip = . (147) 

Another interesting example (especially when discussing inflation) is the scalar field case. There, as we shall see 
in Section 5.4, FI = 0, but in general F 7^ since 8p/ 5p / p/p. Nevertheless, since this case again has only one 
dynamical degree of freedom, we can express the perturbation equations in terms of one single second order equation 
for 'P. In Section 5.4 we shall find the following equation for a perturbed scalar field cosmology 

^ + 'iJf{\+c])'V+[{\+3cl){M'^-K)-{\+3w) ( 2 +K)+k^Yy = 0. (148) 

The only difference between the perfect fluid and scalar field perturbation equation is that the latter is missing the 
factor in front of the oscillatory term. Note also that for K = Q and w = = constant, the time dependent mass 
term m^{7]) = -{\ +T,c]){,ytf^- - K) + {\ +Zw) {M^ + K) vanishes. 

It is useful to define the variable [32] 

M = a[47rG(Jf'2-Jf' + /s:)]"^^V, (149) 

which satisfies the equation 

M+(Y/fe2-0/0)M = O, (150) 
where Y = or Y = 1 for a perfect fluid or a scalar field background respectively, and 

e= . ■ (151) 



Another interesting variable is 



3(1 +w) 
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For the rest of this section we set ^ = for simphcity. Using Eqs. (147) and (148) respectively one then finds 

i; = -e (153) 

On super-horizon scales, k/J^<^ 1, this time derivative is suppressed by a factors {^jM'Y — ik'^f' and this variable 
is (nearly) conserved on large scales. 

The evolution of C, is closely related to the canonical variable v defined by 

v = -^^^^C, (154) 

which satisfies the equation 

v+(Y^^-z/z)v = 0, forz=^^^^^±^. (155) 

5.3. Dust and radiation 

Next we discuss two simple apphcations which are important to understand the anisotropics in the cosmic mi- 
crowave background (CMB). 



53.1. The pure dust fluid for A" = 0, A = 

'Dust' is the cosmological term for non-relativistic particles for which we can neglect the pressure so that w = = 
p = and n = r = 0. The Friedmann equation implies for dust a oc tj^ so tha M' = 2/tj. Equation (147) then reduces 
to 

*+-4' = 0, (156) 

with the general solution 

'P = 'Po + 1'i-!t (157) 

with arbitrary constants 'Pq and ^i. Since the perturbations are supposed to be small initially, they cannot diverge for 
7] ^ 0, and we have therefore to keep only the 'growing' mode, *Pi = 0. But also the 'Pq mode is only constant. This 
fact led Lifshitz who was the first to analyze cosmological perturbations to the conclusions that Unear perturbations do 
not grow in a Friedman universe and cosmic structure cannot have evolved by gravitational instability [33]. However, 
the important point to note here is that, even if the gravitational potential remains constant, matter density fluctuations 
do grow on sub-horizon scales, scales where kr\^\ and hence structure can evolve on scales which are smaller than 
the Hubble scale. 

Defining x = A:?], we obtain for the velocity potential and the density contrast 

y = 'Pol (158) 

Z)„ = -5'Po - i'Pox^ , Z) = D.-h3'P-|--y = -i'Pox^ . (159) 
6 X 6 

In the variable D the constant term has disappeared and we have Z) <c *P on super-horizon scales, x <C 1. 

On sub-horizon scales, the density fluctuations grow Uke the scale factor a «: x^. Nevertheless, Lifshitz' conclu- 
sion [33] that pure gravitational instability cannot be the cause for structure formation has some truth: if we start from 
tiny thermal fluctuations of the order of 10^^^, they can only grow to about 10^^*^ due to this mild, power law instability 
during the matter dominated regime. Or, to put it differently, if we want to form structure by gravitational instability, 
we need initial fluctuations of the order of at least 10"^, much larger than thermal fluctuations. According to what we 
have said here, we need these fluctuations at the beginning of the matter dominated phase, but as we shaU see below, 
perturbations do not grow at aU during the radiation dominated era, so that reaUy initial fluctuations with amplitudes 
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:^ 10~^ are needed. One possibility to create such fluctuations is quantum particle production in the classical gravita- 
tional field during inflation. The rapid expansion of the universe during inflation quickly expands microscopic scales, 
at which quantum fluctuations are important, to cosmological scales where these fluctuations are then "frozen in" as 
classical perturbations in the energy density and the geometry. We wiU discuss the induced spectrum on fluctuations 
in Section 5.5. 



5.3.2. The pure radiation fluid, K = 0, A = 



In this limit we set w = c^ = 1/3, and n = F = so that <J> = —'P. We conclude from p <xa and the Friedmann 



equation that ao< -q. For radiation, the M-equation (150) becomes 



.1 9 2 , 



(160) 



with general solution 



m(x) =A 



sin(x) 



— cos(x) I +B 



cos(x) 



— sin(x) 



where we have set x = kri /\/3 = Cjfcrj . For the Bardeen potential we obtain with (149), up to constant factors. 



(161) 



(162) 



We must set B = for perturbations to remain regular at early times. On super-horizon scales, x <C 1, we then have 

A 



«P(x): 



For the density and velocity perturbations one finds 



D^ = 2A 



cos(x) sin(x) 



In the super-horizon regime, x <C 1, this yields 

A 



D„ = -2A- ^ 
3 * 3\/3 



x\ V = 



2\/3' 



;X . 



(163) 



(164) 



(165) 



On sub-horizon scales, x 1, we obtain oscillating solutions with constant amplitude and with frequency k/ \/3: 



- sin(x) , Dg = 2Acos(x) , I* = — Acos(x)/x^ . 



(166) 



Note that also radiation perturbations outside the Hubble horizon are frozen to first order. Once they enter the horizon 
they start to coUapse, but pressure resists the gravitational force and the radiation fluid fluctuations osciUate at constant 
amplitude. The perturbations of the gravitational potential osciUate and decay Uke 1/a^ inside the horizon. 



5.3.3. Adiabatic initial conditions 

Adiabaticity requires that the perturbations of all contributions to the energy density are initially in thermal 
equilibrium. This fixes the ratio of the density perturbations of different components. There is no entropy flux and thus 
r = 0. Here we consider a mixture of non relativistic matter and radiation. Since the matter and radiation perturbations 
behave in the same way on super-horizon scales, 

D^'^=A+Bx\ Df^=A!+B!x\ y W oc o< x, (167) 
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we may require a constant ratio between matter and radiation perturbations. As we have seen in the previous section, 
inside the horizon (x > 1) radiation perturbations start to oscillate while matter perturbations keep following a power 
law. On sub-horizon scales a constant ratio can thus no longer be maintained. There are two interesting possibilities: 
adiabatic and isocurvature perturbations. Here we concentrate on adiabatic perturbations which seem to dominate the 
observed CMB anisotropics. 

From r = one easily derives that two components with pij pi = wi =constant, i = 1 , 2, are adiabatically coupled if 
{\+w\)Df^ = (1 +W2)I>^'^ Energy conservation then implies that their velocity fields agree, V^^) = V^^). This result 
is also a consequence of the Boltzmann equation in the strong coupling regime. We therefore require 

yW=y('"), (168) 



so that the energy flux in the two fluids is coupled initially. 

We restrict ourselves to a matter dominated background, the situation relevant in the observed universe after equality. 
We first have to determine the radiation perturbations during a matter dominated era. Since is dominated by the 
matter contribution (it is proportional to the background density of a given component), we have ^ ~ const = 'Pq- 
We neglect the contribution from the sub-dominant radiation to Energy momentum conservation for radiation then 
gives, with x = kr], and d/dx = i 

Z)W' = -lyW (169) 

yW' = T¥ + ^D^g\ (170) 

Here ^ is just a constant given by the matter perturbations, and it acts like a constant source term. The general 
solution of this system is then 

D^g^ = Acos(c^jc)- -^Bsin(cjx)-F8'P[cos(c^x)- 1] (171) 
v3 

/3 

yW = £cos(Qx) + ^Asin(c,x)+2V31'sin(c,x), (172) 
where = 1 / \/3 is the sound speed of radiation. Our adiabatic initial conditions require 

y{r) y(m) 

Um = Vo = Um < oo. (173) 

a:->0 X X 

Therefore B = and V = Vqx with Vo=A/4 — 2'P on super horizon scales, x <C 1. Using in addition W = 3Vo (see 
(158)) we obtain 

/ ^ A / Y \ 

(174) 
(175) 

(176) 
(177) 



— 


-•Pcos [ ^ ) - S'P 


yW = 




r,('n) _ 
— 


-^(5 + -x^) 
6 


yW = 





On super-horizon scales, x <C 1 we have 



20 1 
~ -y'P and y W ~ -x'P , (178) 



note that D^g'' = {4/3)Dg"^ and V^''^ = V^""^ as it is required for adiabatic initial conditions. 
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5.4. Scalar field cosmology 



We now consider the special case of a Friedmann universe filled with self interacting scalar field matter. We keep 
spatial curvature ^ = in this section. The action is given by 



S=Y^|d^x^/W\R + |d'xV\g\(^^^,(p^^(p-W{(p)j =Sg+S„ 



(179) 



where (p denotes the scalar field and W is the potential. The energy momentum tensor is obtained by varying the matter 
part of the action, wrt the metric g^^^, 



guv 



The energy density p and the energy flux u axe defined by 

For a homogeneous and isotropic universe, (p = (pit) and g^v = a^r\^y we obtain 



The pressure is given by 



P = ^<i>'+w = ^(1,0) 



We now consider scalar field perturbations. 



(p = (p + 8<p . 



(180) 



(181) 



(182) 



(183) 



(184) 



Clearly, the scalar field only generates scalar-type perturbations (to first order). The perturbed metric is therefore 
given by ds^ = — a^(l + 2^)drf- + a?{\— 2^)dijdx'dx-' . Inserting Eq. (184) in the definition of the energy velocity 
perturbation V, 



and the energy density perturbation Sp, 
we obtain 

and 



p=p + 5p , 

1 ^ 1 ^9 

5p = —(p5<p-^(p'-^ + W,q,5(p 



V = -5(p . 
9 



From the stress tensor, Tij = (p,i (p,j— [^dx (pd^ (p + W] gtj we find 

1 . „ 1 



pnL=—^8(p i^^^m -W,,f,d(p and n = 0. 



Short calculations give 



Dg = -(1+w) 



a 



2W. 



PP 



^ [^pDs-p5(p] 



n = 0. 



(185) 
(186) 
(187) 

(188) 

(189) 

(190) 
(191) 
(192) 
(193) 
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The Einstein equations then lead to the following second order equation for the Bardeen potential which we have 
discussed above: 

'ii' + 2(Jf - 9/9)4'+ (2^-2jr9/(i)+)fc2)'P = (194) 
or, using the definition = p/p, 

* + 3^(l+c?)4'+(2^+(l+3c?),if2 + yt2)^ = . (195) 

As already mentioned above, this equation differs from the ^ equation for a perfect fluid only in the last term 
proportional to k^. This comes from the fact that the scalar field is not in a thermal state with fixed entropy, but it 
is in a fully coherent state (F 7^ 0) and field fluctuations propagate with the speed of Ught. On liirge scales, A:TJ <C 1, 
this difference is not relevant, but on sub-horizon scales it does play a certain role. 

5.5. Generation of perturbations during inflation 

So far we have simply assumed some initial fluctuation amplitude A, without investigating where it came from 
or what the ^-dependence of A might be. In this section we discuss the most common idea about the generation of 
cosmological perturbations, namely their production from quantum vacuum fluctuations during an inflationary phase. 
The treatment here is focused mainly on getting the correct result with as little effort as possible; we ignore several 
subtleties related, e.g. to the transition from quantum fluctuations of the field to classical fluctuations in the energy 
momentum tensor The idea is of course that the source for the metric fluctuations are the expectation values of the 
energy momentum tensor operator of the scalar field. 

The basic idea is simple: A time dependent gravitational field very generically leads to particle production, analo- 
gously to the electron positron production in a classical, time dependent, strong electromagnetic field. 

Let us first fix our notation. Inflation is an era during which the expansion of the scale factor is accelerated, ^ > 0. 
In terms of conformal time, -jfj = , this becomes 

— = -Jf >0 . 
dfi a 

We shall only consider simple power law inflation, where a = (cTj)* for some constants c and q. For the scale factor 
to be positive and real we require cq > 0. Expansion then happens when cq > and accelerated expansion when in 
addition ^ < 0. Hence for power law inflation, the scale factor behaves like 

floc |J7|'? 

and 77 < as well as ^ < 0. It is easy to see that de Sitter inflation, a °= cxp{Ht), corresponds to ^ = — 1. In general, 
for a fluid with p = wp 

2 

1 — 3w 

Inflation therefore requires w < — 1 /3 . During scalar field inflation, the energy density must therefore be dominated by 
the potential, W > a^^<p^. We suppose that the field is 'slowly rolling' down the potential until at some later moment 
the condition w < —1/3 breaks down and inflation stops. How far away a given moment is from this end of inflation 
can be cast in terms of the slow roU parameters fii and £2 defined by 

H Jif 

£1 = 7T , H= — is the Hubble parameter (196) 



£2 



(197) 



5.5.7. Scalar perturbations 

The main result of this subsection is the following: During inflation, the produced particles induce a gravitational 
field with a (nearly) scale invariant spectrum, 

Ar'|'P(A;,Tj)|2 =^-1 x const, with n~ 1 . (198) 
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The quantity k^\^{k,ri)\^ is the squared amplitude of the metric perturbation at comoving scale X=K/k. To ensure 
that this quantity is small over a broad range of scales, so that neither black holes form on small scales nor large 
deviation from homogeneity and isotropy on large scales appear, we must require n ~ 1. These arguments have been 
put forward by Harrison and Zel'dovich [34] (ten years before the advent of inflation), leading to the name 'Harrison- 
Zel'dovich spectrum' for a scale invariant perturbation spectrum. 

To derive the above result we consider a scalar field background dominated by a potential, hence \r]\'^ with 
g ^ — 1 . Developing the action of this system, 

to second order in the perturbations (see [32]) around the Friedmann solution one obtains 

dS = j dx'^\]^{d^vf (199) 
up to a total differential. Here v is the perturbation variable 



(200) 



introduced in Eq. (154). Via the Einstein equations, this variable can also be interpreted as representing the fluctuations 
in the scalar field. Therefore, we quantize v and assume that initially, on small scales, A:]?] | <C 1, v is in the (Minkowski) 
quantum vacuum state of a massless scalar field with mode function 

Vin = -^exp(jA:77) . (201) 

The pre-factor vq is a ^-independent constant which depends on convention, but is of order unity. From (153) we can 
derive 

(v/z) = - , 

where z a is defined in Eq. (155) and u <x. aTj*P is given in Eq. (149). On smaU scales, k\r[ \ <^ 1, this results in the 
initial condition for u 

Min = ^exp(jA:Tj) . (202) 
In the case of power law expansion, a o= \ri\'', the evolution equation for u, Eq. (150), reduces to 

U+{k^-^^i±}l)u = 0. (203) 

The solutions to this equation are of the form {k\r]\y^^H^^^ (k^), where IJ. = q + 1/2 and is the Hankel function 
of the ith kind (i = 1 or 2) of order fj.. The initial condition (202) requires that only i/^^^ appears, so that we obtain 

where again a is a constant of order unity. We define the value of the Hubble parameter during inflation, which is 
nearly constant by Hi. With // = ^/a~l/(|Tj|a), we then obtain a ~ 1 /(//,■ | T] |). With the Planck mass defined by 
4nG = M4 2, Eq. (149) then gives 

On small scales this is a simple osciUating function while on large scales A:| tj | -C 1 it can be approximated by a power 
law, 

'P- ^A:-3/2(^|^|)i+9 ^ for k\r]\ «; 1 . (205) 
M4 
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Here we have used n = l/2 + q <0. This yields 



(^(A:|T?|)i+«) ^^^^^ 

hence n~lif^~— 1. Detailed studies have shown that the amphtude of can still be somewhat affected by the 
transition from inflation to the subsequent radiation era, the spectral index, however, is very stable. Simple deviations 
from de Sitter inflation, like e.g. power law inflation, ^ > — 1, lead to slightly blue spectra, n ~ 1. 

With a somewhat more careful treatment, one finds that both, the amplitude and the spectral index depend on scale 
via the slow roll parameters fii and £2, 

W = T^(^'7/)"-', (207) 
n\k=aMHM = l-4ei(77)-2e2(77). (208) 

Vector perturbations are not generated during standard inflation; and even if they are generated they only decay 
during subsequent evolution and we therefore do not discuss them any further. This may change drastically in 
braneworlds (see [24])! 



5.5.2. Tensor perturbations 

situation is different for tensor perti 

impUes 



(T) 

The situation is different for tensor perturbations. Again we consider the perfect fluid case, 11;^ = 0. Eq. (144) 



Hij + ^Hij + k^Hij = . (209) 

If the background has a power law evolution, a°^r]'> this equation can be solved in terms of Bessel or Hankel functions. 
The less decaying mode solution to Eq. (209) is Hjj = eijX^^'^~I^Ji/2-qix), where Jv denotes the Bessel function of 
order v,x = kr\ and e,y is a transverse traceless polarization tensor. This leads to 

Hij = const for x< 1 (210) 

Hij = - for (211) 
a 

One may also quantize the tensor fluctuations which represent gravitons. Doing this, one obtains (up to log corrections) 
a scale invariant spectrum of tensor fluctuations from inflation: for tensor perturbations the canonical variable is simple 
given by hij = MpaHjj. The evolution equation for hjj = heij is of the form 

h+{k^+m^{r]))h = , (212) 

where = —aja. During inflation = —q{q — l)/??^ is negative, leading to particle creation. Like for scalar 

perturbations, the vacuum initial conditions are given on scales which are inside the horizon, |ot^|, 

h\-a = ^exp(i^T7) for /:|t7| » 1. 



Solving Eq. (212) with this initial condition, gives 



(2) 



(feTj) 



(2) 

where H\, is the Hankel function of degree V of the second kind. On super-horizon scales where we have 
^f-i/i^^^) °^ (^1^1)^"'^^' this leads to ~ |t7|(A:|t7|)2^"^ Using the relation between hij = heij and Hij one obtains 
the spectrum of tensor perturbations generated during inflation. For exponential inflation, ^ ~ — 1 one finds again a 
scale invariant spectrum for Hij on super-horizon scales, 

1^\HijW}\c^{lHi^/Mif{kr]fY^ with nr = 2(^-1- 1) ~ . (213) 
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Again, a more careful treatment within the slow roll approximation gives 

tiT = -2ei . (214) 

A more detailed analysis also of the amplitudes of the scalar and tensor spectra leads to the consistency relation nr = 
~2Aj /A| of slow roU inflation. Here At and As are the amplitudes of tensor and scalar perturbations, respectively. 
More details on inflation can be found in many cosmology books, e.g. Refs. [29, 35, 36]. 



5.6. Power spectra 

The quantities which we have calculated in the previous subsection are not the precise values of e.g. ^{k,ri), but 
only expectation values (|^(k, T])^). In different realizations of the same inflationary model, the 'phases' a(k, ?]) 
given by *P(k, 77) = exp(ia(k))|*I'(/c)| are different. They are random variables. Since the process which generates the 
fluctuations ^ is stochastically homogeneous and isotropic, these phases are uncorrelated (for different values of k). 
However, the quantity which we can calculate for a given model and which then has to be compared with observation 
is the power spectrum. Power spectra are the "harmonic transforms" of the two point correlation functions^. If the 
perturbations of the model under consideration are Gaussian, this is a relatively generic prediction from inflationary 
models, then the two-point functions and therefore the power spectra contain the full statistical information of the 
model. 

Let us first consider the power spectrum of matter, 

7^z)W = (|Z)^(k,T7o)|') . (215) 

Here ( ) indicates a statistical average, ensemble average, over "initial conditions" in a given model. Poik) is usuaUy 
compared with the observed power spectrum of the galaxy distribution. 

The spectrum we can both, measure and calculate to the best accuracy is the CMB anisotropy power spectrum. It 
is defined as follows: The fluctuations of the radiation temperature as observed in the sky, AT/T, is a function of our 
position xo, time tJq and the photon direction n. We develop the n-dependence in terms of spherical harmonics. We will 
suppress the argument r/o and often also xq in the following calculations. All results are for today (tjq) and here (xq). 
By statistical homogeneity statistical averages over an ensemble of realizations (expectation values) are supposed to 
be independent of position. Furthermore, we assume that the process generating the initial perturbations is statisticaUy 
isotropic. Then, the off-diagonal correlators of the expansion coefficients agm vanish and we have 

^ (xo, T}o,n) = J^a£;„(xo, 7]o)I'^m(n), {aem ■ a\,J) = Su'S^'Q ■ (216) 

The Q's are the CMB power spectrum. 

The two point correlation function is related to the Q's by 

^(n)^(n')\ = £ (a,^.4,,)y,^(n)y;„,(n') = 

/ n-n'=^i e,e',m,m' 
^ 1 

£ y^^(n)y4(n') = — £(2^+ 1)QPKm), (217) 

e m=-e e 



where we have used the addition theorem of spherical harmonics for the last equaUty; the Pi's are the Legendre 
polynomials. 

For given metric perturbations and perturbations of the energy momentum tensor of the cosmic fluid, the tem- 
perature perturbations can be determined by following the oscillations in the radiation fluid before decoupling (see 



^ The "harmonic transform" in usual flat space is simply the Fourier transform. In curved space it is the expansion in terms of eigenfunctions of the 
Laplacian on that space, e.g. on the sphere it corresponds to the expansion in terms of spherical harmonics 
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subsection 5.3.2) and by following the propagation of photons along geodesies in the perturbed spacetime after de- 
coupling. Decoupling of photons and matter happens during recombination (T ~ 3000K, z — 1000), where electrons 
and protons recombine to neutral hydrogen. During that process, the number density of free electrons with which the 
photons can scatter drops drastically and finally becomes so low, that the mean free path of the photons grows larger 
than the Hubble scale. The surface of constant temperature, T^cc = 7"('7dcc)> ^it which this happens is also called the 
'last scattering surface' . After last scattering, the photons effectively cease to interact and move freely along geodesies 
(more details can be found e.g. in [29, 28, 37]). 
Clearly the a/^'s from scalar-, vector- and tensor-type perturbations are uncorrelated. 



Since vector perturbations decay, their contributions, the Q % are negligible in models where initial perturbations 
have been laid down very early, e.g. , after an inflationary period. Tensor perturbations are constant on super-horizon 
scales and perform damped oscillations once they enter the horizon. 

Let us first discuss in somewhat more detail scalar perturbations. We restrict ourselves to the case K = for 
simplicity. We suppose the initial perturbations to be given by a spectrum. 



We multiply by the constant rj^^ , the actual comoving size of the horizon, in order to keep A dimensionless for all 
values of n. A then represents the amplitude of metric perturbations at horizon scale today, k= I/170. 
For adiabatic perturbations we have obtained on super-horizon scales. 



The dominant contribution to the temperature fluctuations on super-horizon scales (neglecting the integrated Sachs- 
Wolfe effect /<i> — 'i' ) comes from two terms: the first, 2^, is the change of photon energy due to the gravitational 

(r) 

potential at the last scattering surface, 7] = 77dec, and the second, ^/^Dg represents the intrinsic temperature fluctuations 
(for more details see [38, 30, 37]). With Eq. (220) this yields the famous Sachs- Wolfe formula 




(218) 




(219) 




(220) 



— (xo,n,T]o) =2'P(xdec,I7dec) + 4£*r(-«dec,I7dec) = ^ '^'(•^dec) ^Tdec)- 



(221) 



The Fourier transform of (221) gives 



(222) 



Using the decomposition 



^,k„(^o-')dec) = ^ (2£ + \)fk{k{r]o - T?dec))P^(k- n) 



where are the spherical Bessel functions, we obtain 




(223) 



■iiiKno - nito))je (Kno - ■ndec))Pe{^ ■ n) • P^'ik ■ n') . 



(224) 
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In the second equal sign we have used the unitarity of the Fourier transformation. Inserting P({kn) = 
^'LmYi!',„ik)Yem{n) and Pfikn') = 2F+Tl^'n'^fm'(^)^^''n'("'). the integration over directions dQ.^ gives 

5«'5„„„'i:,„}'4(n)}'fo,(n'). 

Using as well 'LmYemi'^)^i'ni'a') = ^Peil^^), where ;U = n • n', we find 

^(xo,n,?7o)^(xo,n',T7o) 

^ ^ ' nn'=M 



(^I'l'rt^^yfW^o-W)- (225) 



Comparing this equation with Eq. (217) we obtain for adiabatic perturbations on scales 2 < ^ -C (t]o — TJdec)/T?dec 
■100 

'■ 1 



n Jq k 



3^ 



k'jfikino-ndec)). (226) 



If 1* is a pure power law as in Eq. (219) and we set k{rio — rjdeo) ~ krio, the integral (226) can be performed 
analytically. For the ansatz (219) one finds 

Asw) r(3-»)r(£-i + §) n-yi, 

C) ' = — — — — - for — 3 < n < 3 . (227) 

i 9 23-«r2(2-^)r(^+§-f) 

Of special interest is the scale invariant or Harrison-Zel'dovich spectrum, n = 1 (see Section 5.5). It leads to 

^(£+l)c(^)=const.~^(^^(t?^)j ^ , ^i = K/t. (228) 

This is precisely (within the accuracy of the experiment) the behavior observed by the DMR experiment aboard the 
satellite COBE [39] and more recently with the WMAP satellite [40]. 

Infiationary models predict very generically a HZ spectrum (up to small corrections). The DMR discovery has 
therefore been regarded as a great success, if not a proof, of inflation. There are other models like topological defects 
[41] or certain string cosmology models [42] which also predict scale-invariant, i.e. Harrison Zel'dovich spectra of 
fluctuations. These models do however not belong to the class investigated here, since in these models perturbations 
are induced by seeds which evolve non-Unearly in time. 

For gravitational waves (tensor fluctuations), a formula analogous to (227) can be derived. 



- [ dkk^l driH(ri,k) 



^ ^ . (229) 



To a very crude approximation we may assume H = 0on super-horizon scales and / drjH ji{k{r}Q — J])) ~ H{ri = 
l/k) je {krjo). For a pure power law. 



one obtains 



kU\Hik,r] = l/k)\') ^A'^k"^r^,"\ (230) 



C<J) ^ 2 (£ + 2)! ^2 fdx^rjM 

TJ ^ 



n{£-2)\ 'J X x^ 
{£ + 2)1 . r(6-nr)r(f-2 + f ) 



{£-2)1 ^26-«rr2(Z-„y)r(^ + 4-^^)' 
For a scale invariant spectrum (nr = 0) this results in 



(231) 



^(D _ £{£+!) 

^ {£ + 3){£-2)"' 15;r 



• + 1)Q ' ^ , .\T; .A t^: ■ (232) 
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The singularity at ^ = 2 in this crude approximation is not real, but there is some enhancement of + 1)Q at 
£ ^ 2 see Fig. 7). Again, inflationary models (and topological defects) predict a scale invariant spectrum of tensor 
fluctuations (ut ~ 0). 

On intermediate scales, 100 < f < 1000, the acoustic oscillations of radiation density fluctuations before decoupling 
(see subsection 5.3.2) lead to a characteristic series of peaks in the CMB power spectrum which is being measured 
in great detail and contains very important information on cosmological parameters [43]. On small angular scales, 
£ ~ 800, fluctuations are damped by coUisional damping (Silk damping). This effect has to be discussed with the 
Boltzmann equation for photons, which goes beyond the scope of this introduction (see [29, 43]). 



BxlO-'o 

Bxio-'" 
txio-'» 

SxlO-'o 



10 100 1000 10 100 1000 





10 100 1000 10 100 1000 

1 1 

FIGURE 7. Adiabatic scalar (left) and tensor (right) CMB anisotropy spcetra are shown. The dimensionless quantity £{£ + 
l)Q/(2;r) is plotted. The top panels show the temerature anisotropies while the bottom panels show the corresponding polarization 
spectra (for an introduction to polarization see e.g. [29]). 



6. BRANEWORLD COSMOLOGY 

We now want to study cosmology of an expanding maximaUy symmetric braneworld. We still require the bulk to be 
empty and Z2-symmetric. One can show that the most general empty bulk allowing for a homogeneous and isotropic 
brane is Schwarzschild-AdS (Sch-AdS). In the cosmological setting we allow the brane to move in the bulk. As we 
shall see, this can mimic cosmological expansion. The situation is as depicted in Fig 8. 
The 5d metric of Sch-AdS is of the form 

ds' = -F{R)dT' + 1^ + ^'^^') (233) 

where the function F is determined by the AdS curvature radius £, the 5d mass C and the curvature K of 3d space (on 
the brane) via 
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FIGURE 8. A moving brane in AdS-Schwarzschild. In this coordinates the AdS horizon is at /f = 



From this we can calculate the 5d Weyl tensor, Cabcd, and its 'electric' components defined in (30) with the result 

E^y = p'^u^u^ + K^y with M = 1, 0) , (234) 

Eoo = ^=P^ ?rfv = 0. (235) 

Here R{T) = a{t), where t is cosmic time on the brane and u is the unit normal vector in direction of cosmic time on 
the brane. R{T) is the brane position at time t{T). 



6.1. The modified Friedmann equations 

From the brane gravity equations, 

with A4 = 2 (A5 + JCs /6X^), K4 = K-|A /6 and 
with T^v = (P + p)UfiUv - pg^v we obtain 

where H is the Hubble parameter. In this section we denote the derivative with respect to cosmological time t 
determined by dt = adrj by an over-dot, so that H = d/a and Jf = d. 

The term p^/ (2A) in Eq. (236) is a correction to the Friedmann equation which is important only at high energies 
and = Cja^, comes from the Weyl tensor. It is called 'Weyl radiation' since it scales like cosmic radiation, 
oc a~'^. Observations (nucleosynthesis) tell us that latest at the temperature T IMeV these additional terms should 
be unimportant. More precisely, P^/p| ^ 0-1 and A ~ (IMeV)"^. For the 5d Planck mass this implies M5 ~ 
3 X lO^GeV. The conservation equation of 4-dimensional cosmology remains unchanged, 

P = -3(p + p)-. (237) 
a 

Solutions for C = = A4 = are readily found. If the equation of state is given by p = wp we find 

a = ao[/(/ + ?A)]^ , tx= -<lsec, (238) 

v3;rA(3 -l-w) 
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where we have used A > l(Mev)'* for the inequahty. This is to be compared with the usual 4-dimensional behavior 
from general relativity (GR). There we have a = aot^^^^^^"\ which corresponds to the above result in the limit X 



Especially interesting is also the case of de Sitter inflation where p ■ 
with 



-p and hence p =constant, so that a = a^e 



Ht 



H- 



> 



— =Hgr 



(239) 



At low energies we recover the usual Friedmann equation while at high energies, p > A, the expansion law differs. 
During a radiation epoch, p = p/3 with p » A we have a <x t^/^ (instead of the usual GR behavior, a x t^l'^). This 
comes from the fact that in braneworlds at high energies H <^ p, where as in 4d GR we have H <^ If A is given 
by the electroweak scale, A > ITeV*, the observed low energy cosmology, hke nucleosynthesis which starts after 
p ~ IMeV* is not affected (if C is sufficiently small). 

However, perturbations will carry 5d effects which should in principle be observable in the fluctuation spectrum of 
the cosmic microwave background radiation (CMB) and in the large scale distribution of matter. These effects are still 
under investigation. In Section 6.3 I shall present some prehminary partial results. 



6.2. Brane inflation 

We now want to study scalar field inflation in braneworlds. Since energy momentum conservation is still valid, the 
scalar field evolution equation is also not modified, 

(p + 2>H(p + W'{(f>)=Q . 

In 4d GR, the condition for inflation, d > , is equivalent to 9^ < V which corresponds to the violation of the strong 
energy condition 



1 



W<-\p 



1 



< 



1 

3 ' 



or q)^ <W . 



The braneworld Friedmann equation (239) leads to a stronger condition on w for accelerated expansion. For 
branewords < ^=H + requires 



1 

w<-- 



l+2p/A 



1+p/A 



or q)'^ <W- 



A 



4*^ 



1 

2^ 



(240) 



for inflation to happen. This becomes the usual w < — 1/3 at low energy, p << A, but turns into w < —2/3 at high 
energy. 

If the slow roll approximation (9^ <c IV) is satisfied we have 



K4 



W 



1 + 



w 

2A 



'3H 



Since the Hubble rate is increased, slow roll is maintained longer than in usual 4d inflation. Correspondingly, the slow 
roll parameters are reduced. 



_ H _ Ml 





\+w/x 







Mi 
8;r 



(^) 


1 




X\+W/2X)_ 



(241) 



In the high energy regime, V » A they are reduced by factors 4A /V and 2A /V, respectively. Hence, the universe may 
be inflating only because it is in the high energy regime and turn into kinetic dominated expansion, w ~ 1, as soon as 
y < A. Such models are constrained since they induce a blue spectrum of gravity waves [44]. 

In standard 4d GR the perturbation spectrum induced by inflation is well known and the scalar spectrum agrees 
extremely well with the observed anisotropies in the CMB. This will most probably lead to the most stringent 
constraints for braneworlds, which however have not yet been explored in fuU generality so far. This is still a very 
active field of research. 
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6.3. Observable consequences from braneworld cosmology 



So far, we have seen that it is conceivable that our Universe is a 3-brane. At least at low energy and disregarding 
perturbations, we cannot distinguish cosmological evolution due to 4d Einstein equations or the (so different!) brane 
gravity equations. We finally want to study ways to discover whether the braneworld idea is realized in nature. Are 
there tell-tale observational signatures which would betray whether we live on a brane? 

As we have seen, the gravitational equations for braneworlds differ significantiy from Einstein gravity. However, at 
low energy the Friedmann equations for a homogeneous and isotropic Universe are recovered. Hence there are two 
regimes in which deviations from Einstein gravity will be found: 

• At high energy. This is especially interesting for the generation of inflationary perturbations which may 
be affected by high energy braneworld behavior. As long as we restrict ourselves on background effects the 
calculations are relatively straight forward and well under control. 

• In the perturbations. Cosmological brane perturbation theory is not well under control and still a subject of 
active research. One main point iire bulk perturbations which can affect the brane and act there like 'sources'. On 
the other hand, gravity wave perturbations generated on the brane can be emitted into the bulk. 

Here we give examples of both aspects, how effects from braneworlds can enter cosmological perturbations, but we 
are by no means exhaustive (more details and especially references can be found in [16]). 

Let us first consider the high energy universe. During inflation scalar and tensor perturbations are generated. The 
spectrum of scalar perturbations is I'PpA:^ = Agk.^i~^ = A|(A:///o)"* The slow roll approximation for braneworlds 
gives [16] 



o 5l2n 
As 



75M| W'^ 



ns = 1 - 4ei - 2e2 . (242) 



2A 

Similarly, for tensor perturbations, \H\^k^ = Ajk^ one obtains in the slow roll approximation [45] 
SW_ 



A\ ~ ^^^F^{H/iJ.) , nr = -2ei, where F{x)= Vl+x^-x^smh~\l/x) 




-1/2 

(243) 

Combining these slow roll equations, at low energy one obtains the same consistency relation as for ordinary inflation. 

At higher energies, however the relation is different. Furthermore, the tensor to scalar ratio R = (At/As)^ and the 
spectral indices, both depend on the energy scale W. In Fig. 9 the behavior of the different quantities is indicated as 
function of the energy. Of course also the amplitude of the perturbations strongly depends on the parameter W^/M4. 

In Fig. 10 two models are shown, quartic inflation with W = acp'^ and quadratic inflation with W = rrp'(p^. The lines 
of these models in the {R,ns) plane for varying A are drawn. The parapeters m respectively a are chosen such that 
the scalar amplitude is compatible with the measured value, A| ~ 10"^*' for each brane tension A. The observational 
constraints from WMAP data [40] are also indicated. It is clear, that quartic braneworld inflation fares even worse than 
ordinary quartic inflation. It is virtually excluded. Also quadratic inflation with strong braneworld effects fits the data 
somewhat less well than usual quadratic inflation, since it predicts too strong tensor contributions. But clearly, in lack 
of a concrete model of inflation (e.g., a given potential) there is Uttle which can be said. 

Discussing the effects on perturbations from braneworlds is opening Pandora's box. There is a plethora of new 
phenomena some of which we don't even know the sign. For example: during ordinary inflation, gravitational 
waves are generated. For a given inflationary potential, their amplitude can be calculated accurately. However, in 
the braneworld context, a fraction of these waves will be radiated into the bulk and thereby reduce the gravity wave 
amplitude. On the other hand, there is also gravity wave generation in the bulk, and some of these accumulate on the 
brane, increasing the amplitude of gravity waves on the brane. Therefore, depending on the precise realization, even 
the sign of the braneworld effect on a gravity wave background is unknown. 

For a more concrete example, let us concentrate on scalar perturbations. We just take into account, that on the 
perturbative level the Weyl tensor E^y can no longer be neglected. Its energy density perturbation 5pE acts Uke 
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FIGURE 9. We show the slow roll parameters £i (dotted, red) and £2 (long dashed, blue) as function of the energy scale of 

inflation, V /X . The spectral index n (solid, black) and the tensor to scalar ratio R (dashed, magenta) are also indicated. For IV/A 1 
the slow roll parameters stay nearly constant and correspond to their initial values which are chosen £1 (0) = 0. 1 and £2(0) = 0.02. 



Ouartic \ \ 

h.e.-^A</l.e 


/ 


! [ . i A 


i 



-0.1 -0.09 -0.08 -0.07 -0.06 -0.05 -0.04 -0.03 -0.02 -0.01 
n-1 



FIGURE 10. The dependence on the brane tension X of the parameters {R,ns) is shown for quartic and quadratic inflation. The 
1-, 2- and 3-a contours from the WMAP experiment are also indicated (from [45]). 



a radiation perturbation. In addition, however it can have an arbitrary amount of anisotropic stress. He- The latter 
induces a difference between the two Bardeen potentials, ^ — 4) oc n^. This affects mainly the Sachs- Wolfe term in 
the CMB fluctuations, hence the low multi-poles up to roughly the first peak. In Fig. 1 1 we show the effect of a Weyl 
perturbation as function of an amplitude parameter 

Cdark = —7J ■ (246) 

Here C,m is the ^ variable defined in Eq. (155), due to ordinary matter (without the Weyl component). The anisotropic 
stress , on large scales, £ ^ Xjk, can be determined as function of Cdark and ^„,. Confidence plots for the amplitude 
Cdark and several other cosmological parameters from the WMAP data are shown in Fig. 12. 

There are many more effects which may come from perturbations in the bulk and the different perturbation equations 
on the brane which are presently under study. A systematic investigation is still lacking. 
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FIGURE 11. The dependence of the CMB power spectrum on the amplitude of the Weyl perturbation, Qark for a fixed set of 
other cosmological parameters corresponding to the concordance model (from [46]). 
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FIGURE 12. Confidence plots from the CMB data for the amplitude Qajk ™d other cosmological parameters. The line contours 
indicate 2- and 3-a. Clearly, a Weyl contribution to the perturbations of more than about 10% is strongly disfavored by the data. 



7. CONCLUSIONS 

In these lectures we have studied the possibility that our Universe may represent a 3-brane in a higher dimensional 
space. This idea is motivated by string theory. We have especially investigated the case of one large extra dimension 
where the brane gravitational equations can be obtained purely from the bulk equations. Even though the resulting 
gravity on the brane differs strongly from Einstein gravity, we have seen that for an Anti-de Sitter bulk, Newton's 
law is recovered at large distances and the Friedmann equations for the evolution of the Universe are obtained at low 
energy. 

It is, however by no means clear, to which extent the different gravitational equations will spoil the successes of 
cosmological perturbation theory. This is still an open question and its answer will be crucial for braneworlds. 
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